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ABSTRACT 

Using the South Pole Telescope (SPT), we have discovered the most massive known galaxy cluster 
at 2; > 1, SPT-CL J2106-5844. In addition to producing a strong Sunyaev-Zel'dovich effect signal, 
this system is a luminous X-ray source and its numerous constituent galaxies display spatial and color 
clustering, all indicating the presence of a massive galaxy cluster. VLT and Magellan spectroscopy 
of 18 member galaxies shows that the cluster is at z = 1.132^g'QQ3. Chandra observations obtained 
through a combined HRC-ACIS GTO program reveal an X-ray spectrum with an Fe K line redshifted 
by 2; = 1.18 ± 0.03. These redshifts are consistent with the galaxy colors found in extensive optical, 
near-infrared, and mid-infrared imaging. SPT-CL J2106-5844 displays extreme X-ray properties for 
a cluster, having a core-excluded temperature of kT = ll.Ol^^g keV and a luminosity (within rsoo) of 
Lx(0.5 — 2.0 keV) = (13.9 ± 1.0) x 10^* erg s~^. The combined mass estimate from measurements of 
the Sunyaev-Zel'dovich effect and X-ray data is M2m = (1.27 ±0.21) x lO^''^ /j.^ Mq. The discovery 
of such a massive gravitationally collapsed system at high redshift provides an interesting laboratory 
for galaxy formation and evolution, and is a probe of extreme perturbations of the primordial matter 
density field. We discuss the latter, determining that, under the assumption of ACDM cosmology 
with only Gaussian perturbations, there is only a 7% chance of finding a galaxy cluster similar to 
SPT-CL J2106-5844 in the 2500 deg^ SPT survey region, and that only one such galaxy cluster is 
expected in the entire sky. 

Subject headings: galaxies: clusters: individual (SPT-CL J2106-5844) galaxies: formation — galax- 
ies: evolution — early universe — large-scale structure of universe 
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1. INTRODUCTION 

Galaxy clusters are the most massive collapsed ob- 
jects in the Universe. Massive clusters form from the 
most overdense perturbations on their mass scale in the 
primordial matter density field and are exceedingly rare 
objects. At high redshift, corresponding to a relatively 
short time after the Big Bang, they provide a power- 
ful probe of the initial matter density field of the Uni- 
verse. Their abundance is a particularly sensitive probe 
of Gaussianity in the rn atter density field (e.g., Fry 1986; 
Mortonson et al.|[2010 |, and, as a result, cluster surveys 
can test va rious models of infl ation and topological de- 
fects (e.g., Verde et al. 2001). Observations of galaxy 



clusters, in which the constituent galaxies have similar 
formation histories, ca n strongly constrain hiera rchical 
formation models (e.g., De Lucia fc Blaizot] 2007); high- 
mass clusters at high redshift are particularly interesting 
astronomical laboratories for galaxy formation and evo- 
lution. 

The bulk of baryonic matter in galaxy clusters is in 
the form of fully ionized intracluster gas. The free elec- 
trons in this gas emit thermal Bremsstrahlung radiation 
in the X-ray, and lead to inverse Compton scattering of 
the cosmic microwave background (CMB) photons. This 
scattering results in a spectral distortion of the CMB ob- 
served toward the clus ter know n as the thermal Sunyaev- 
Zel'dovich effect (SZ; [Sunyaev fc' Zerdovich||1972l ) . The 
surface brightness of the ti'A effect is determined solely 
by the thermal pressure of the intra-cluster medium in- 
tegrated along the line of sight and is independent of the 
cluster redshift. The total SZ signal of a cluster is a mea- 
sure of its total thermal energy and is therefore tightly 
correlated with the cluster mass. This makes the SZ 
effect a powerful technique for discovering galaxy clus- 
ters at high redshift and an efficient means of find ing the 
highest-mass clusters (e.g., Carlstrom et al.|[2002| . 

With the goal of detecting galaxy clusters using 
SZ effect, the 10 -m South Pole Telescope (SPT; [Carl 



the 



strom et "aL||2009 ) is surveying 2500 deg^ in three bands 
at millimeter wavelengths. The SPT survey is sensitive 
to galaxy clusters with a mass > 5 x 10^^ Mq at all 
redshifts. The first three SZ-discovered clusters were 



repor ted by the SPT collaboration ( Staniszewski et al 
2009| . The SPT coUaboration recently pubhshed the tuff 
5-(T cluster catalog from the first 200 deg^ o f the survey, 



observed du ring the 2008 observing season (Vanderlinde 
et al. 2010). Most recently, a sample of the 26 most 
massive galaxy clusters extracted from the full 2500 deg^ 
SPT survey has been used to test th e standard cosmolog- 
ical model ( Williamson et aL]|2011 |. We refer the reader 
to these previous SP'l' papers tor a complete description 
of the survey strategy and goals. The Atacama Cos- 
mology Telescope (ACT) has also recentl y reported their 
initia l catalog of SZ-discovered clusters ( Marriage et al. 
20T0|. 



The discovery of z > 1 galaxy clusters has been chal- 
lenging. Recently, several groups have detected galaxy 
clusters at z > 1.1, with most either being found through 
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X-ray 



;e.g. 



Rosati et al.l2004||Mullis et al.|2005 



et al.||2006| |Rosati et al. 20 D9|) or Spitzer mid-intrared 

:-d et 



Stanford 



TMiR) imaging (e.g., Stanfor d et al.|20 "05^, 'Brodwin et al 
2006; Eisenhardt et al. 2008; IMuzzin et al. 2009; WilsoS 
et al.||2 009 Papovich et al. 2010'). Only one SZ-detected 



discovered galaxy cluster has been spectroscopicall y con 
firmed at z > 1 prior to the results presented here ( Brod 



win et aL]|2010 ). In principle, high-mass clusters at high 



redshift can provide evidence for an incomplete under 
standing of cosmology, much like the massive clusters at 
z > 0.3 discovered in the late 1990s provided evidence 
against the theoreticall y preferred = 1 cosmological 



model of th e time (e. g., | Carlberg et a l.|1996 Henry|1997 



BahcaU & Fan 1998 Donahue et al. mm 



During the 2009 campaign of the SP'l-S/ survey, SPT- 
CL J2106-5844 was detected with a signal-to-noise ratio 
(S/N) of 18.5 in 150 GHz data alone (S/N = 22.1 in 
multi-band data) — the largest S /N of any cluster in the 
2008-2009 SPT-SZ sample. Imtial deep optical follow-up 
observations showed no obvious clustered galaxies in gri, 
but did show a strong detection of galaxies in the z band. 
Additional near-infrared and Spitzer photometry further 
confirmed the presence of clustered galaxies. These ex- 
tremely red galaxies are consistent with a galaxy cluster 
at z > 1. 

Chandra observations of SPT-CL J2106-5844 revealed 
a luminous and extended X-ray source. Optical spec- 
troscopy of member galaxies shows that the cluster is at 
z = 1.132. SPT-CL J2106-5844 is therefore not only one 
of the most massive SPT clusters, it is also the highest- 
redshift spectroscopically confirmed SPT cluster. In fact, 
it is the most massive cluster known at redshift z > 1, 
making SPT-CL J2106-5844 an extremely interesting ob- 
ject for detailed study. 

We present our initial detection and follow-up observa- 
tions of SPT-CL J2106-5844 in Section [2| In Section [3) 
we show that SPT-CL J2106-5844 is a > 10^^ Mq galaxy 
cluster at z = 1.132. We examine the implications of the 
existence of such a massive cluster at such high redshift 
m Section H We summarize and conclude in Section [H 
Throughout this paper, Af2oo and Afsoo masses are de- 
fined as the mass enclosed in a spherical region which 
has a density 200 and 500 times the mean and critical 
density of the Universe, respectively. Uncertainties are 
reported for 68% confidence intervals. 

2. OBSERVATIONS, DATA REDUCTION, & INITIAL 
FINDINGS 

2.1. Microwave Observations from the South Pole 
Telescope 

Currently the SPT has been used to complete observa- 
tions of 1500 deg^ to full survey depth of 18 /zK-arcmin 
at 150 GHz and and additional 1000 deg^ has been sur- 
veyed to a "preview" depth of 54 /iK-arcmin, where the 
unit K refers to equivalent fiuctuations in the CMB tem- 
perature. By the end of the 2011 observing season, we 
expect the full 2500 deg^ survey region to be complete 
to full depth. Details of th e survey strategy an d data 
processing are described in iStaniszewski et al.l (f2009), 
Vanderhnde et al.| ( |20To| , and Wilfiamson ct al. (2011). 
I'he 2009 observations that yielded the discovery of SPT- 
CL J2106-5844 included measurements at 95, 150, and 
220 GHz to a noise level of 42, 18, and 85 /xK-arcmin. 

SPT-CL J2106-5844 was detected in maps of 2009 
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SPT data created using time-ordered data processing and 
map-making procedures e quivalent to those described in 
Vanderlinde et al. (2010). Two different cluster extrac- 
tion pipelines were used to detect and characterize clus- 
ters in these SPT maps. Both pipelines use a multi-scale 
matched-filter app roach similar to the one described in 
Melin et al. ( 2006 1 . One pipeline, used to extract clusters 
from single-band 150 GHz data, is id entical to the clus- 
ter ex traction procedure described in Vanderlinde et al. 
(2010). The other pipeline, used to extract clusters from 
multi-band SPT data is identical to that described in 



Williamson et al. (2011). The final SZ observable pro- 
duced by both pipelines is the detection significance 
maximized across all filter scales. SPT-CL J2106-5844 
was detected at ^ = 18.5 by the single-band pipeline and 
at ^ = 22.1 by the multi-band pipeline. 

An image of the filtered, multi-band SPT data — us- 
ing the filter that maximized ^, and centered on SPT- 
CL J2106-5844 — is shown in the left panel of Figure [ll 
Detection significance contours from this same filtered 
map are overlaid in Figures [l] and [1] 

2.2. Optical and Infrared Imaging 

In June 2009, we obtained griz imaging of SPT- 
CL J2106-5844 with the Inamori M agellan Areal C am- 



era and Spectrograph (IMACS; Dressier et al. 



Low Dispersion Survey Spectrograph (LDSSJ ) on the 
twin Magellan 6.5 m telescopes . We obtained additional 



2006 ) and 



/-band imaging with F0RS2 ( [Appenzeller et al.|[T998]) 



on the VLT 8 m telescope in September 2010. See High 
|et al. (2010) for a description of our observation strategy 
and reduction procedure. 
SPT-CL J2106-5844 was observed at the CTIO 4 m 



Blanco tel escope using the NEWFIRM imager (Autry 



et al.|2003 1 on 28 July 2010. Data were obtained in the J 
and Kg filters under photometric conditions. All images 
were acquired using a 30 point random dither pattern. 
At each dither position, six frames with 10 s exposure 
times were coadded, providing 60 s total exposure per 
position. The dither pattern was repeated as necessary 
to achieve total exposure times of 7200 s in J and 4440 s 
in Ks- 

NEWFIRM data were reduced using the FATBOY 
pipeline, originally developed for the FLAMINGOS-2 in- 
strument, and modified to work with NEWFIRM data in 
support of the Infrared Bootes Imaging Survey (Gonza- 
lez, in prep.). Individual processed frames are combined 
using SCAMP and SWARP (Bertin et al. 2002 ), and pho - 
tometry is calibrated to 2MASS (Skrutskie et al. 2006). 
The final images in both filters have FWHM of 1 

Mid-infrared Spitzer/IRAC imaging was obtained in 
September 2009 as part of a larger program to follow up 
clusters identified in the SPT survey. The on-target ob- 
servations consisted of 8 x 100 s and 6 x 30 s dithered 
exposures at 3.6 and 4.5 ^m, respectivelyp^ The deep 
3.6 /im observations are sensitive to passively evolving 
cluster galaxies down to O.IL* at z — 1.5. The data 



were reduced following the method of |Ashby et alT] (2009 ) 
Briefiy, we correct for column pulldown, mosaic the indi 



^'^ http:/ /www.lco.cl/telcscopcs-information/magcUan/ 
instruments-l/ldss-3-1/ . 

In this paper, we will refer to photometry in these bands as 



vidual exposures, resample to 0"86 pixels (half the solid 
angle of the native IRAC pixels), and reject cosmic rays. 

Images generated from the OIR data are presented in 
Figure T] From the deep optical-only images (gri), it is 
difficult to see any indication of a cluster. However, the 
addition of infrared data shows a clear over-abundance of 
very red galaxies clustered near the peak of the SZ decre- 
ment. The over-abundance of galaxies that are spatially 
coincident and have similar colors is a clear indication of 
a cluster of galaxies. Color-magnitude diagrams (CMDs) 
for galaxies within the imaging field of view (FOV; Fig- 
ure p. ) show an overabundance of very red galaxies and 
a probable red sequence. These data are consistent with 
the spectroscopic redshifts described below. 

2.3. Optical Spectroscopy 

We obtained low-resolution spectra of several clus- 
ter members of SPT-CL J2106-5844 with the Gladders 
Image-Slicing Muhislit Option (GISMO; Gladders et al., 
in prep.) module on IMACS mounted on Magellan 
in June 2010; however, these spectra did not yield a 
definitive cluster redshift. Magellan slit masks were de- 
signed using galaxies selected from deep optical and MIR 
imaging compatible with the red sequence expected at 
z > 1.2, as this was the initial redshift approximation 
fr om the photometry. Th e selection was similar to that 
of Brodwin et al. (20101. This technique gives higher 
priority to brighter galaxies with colors consistent with 
the red sequence. If there is a region without a suit- 
able red-sequence target, blue galaxies were selected with 
the idea that they may be emission-line cluster galax- 
ies. The Magellan observations were obtained using the 
f/2 camera with the the 300 1/mm "red" grism and the 
WB6300-9500 fiher. A total of sixteen 1800 s exposures 
were obtained on June 6 and 7. 

The VLT observations were obtained with the 
GRIS_300I grism combined with the OG590 blocking fil- 
ter. Two VLT/F0RS2 masks were designed using galax- 
ies selected from deep / band VLT pre-imaging and the 
Spitzer IRAC imaging to have colors compatible with the 
red sequence expected at z 1.2 to 1.3, which was ini- 
tially indicated by the photometry. Each mask was pop- 
ulated with approximately 30 slitlets of typical length 
10". This target was originally accepted as Priority B 
for VLT/F0RS2 spectroscopy as part of a larger SPT 
spectroscopy campaign. In late November 2010, Chan- 
dra X-ray spectroscop y ga ve an initial redshift from the 
Fe K line (see Section 2.4). After this result when it be- 



[3.6] and [4.5], respectively. 



came clear that there was a significant chance of it not 
being observed in Fall 2010, we submitted a special re- 
quest to the ESO Director to raise this target to Priority 
A. That proposal was successful, and on each of three 
nights in early December 2010 a single 2700 s exposure 
was obtained for a total exposure time of 8100 s on one of 
the two masks. These observations were taken under ex- 
cellent seeing conditions (0'.'5 to 0"8) but at high airmass 
between 1.5 and 1.7. These data were immediately re- 
duced to reveal redshifts of 15 cluster members, enabling 
a robust cluster redshift. 

Spectroscopic reductions proceeded using standard 
CCD pro cessing with I RAF and the COSMOS reduction 
package (Kelson 2003) for the VLT and Magellan data, 
respectively^ 'I' he data were extracted using the opti- 



mal algorithm of Horne ( ,1986) . Flux calibration and tel- 
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Fig. 1. — SPT-CL J2106-5844 at millimeter, optical, and infrared wavelengths. Left: The filtered SZ significance map derived from 
multi-band SPT data. The frame subtends 12' X 12'. The negative trough surrounding the cluster is a result of the filtering of the time 
ordered data and maps. Right: LDSS3 optical and Spitzer/IRAC mid-infrared gi[3.6] (corresponding to BGR channels) images. The frame 
subtends 4f8 X 4.'8. The white contours correspond to the SZ significance from the left-hand panel. The circles mark spectroscopically 
confirmed cluster members, where green indicates quiescent, absorption-line member galaxies and cyan indicates an active, emission-line 
member galaxy. Some spectroscopic member galaxies are outside the FOV for this image. 
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Fig. 2. — Color-magnitude diagram (J — [3.6] vs. [3.6]) for galax- 
ies within the IR AC FOY. Suspected red-sequence cluster members 
are plotted in red. Lower-probability, but potential cluster mem- 
bers are plotted in blue. Spectroscopic members are plotted as 
stars, where the red stars correspond to passive galaxies and the 
blue star represents an emission-line galaxy. Additional galaxies 
in the field are plotted as black points. The size of the symbol is 
inversely proportional to the distance to the center of the cluster 
as determined by the clustering of the red-sequence galaxies. Our 
5-a- limits are plotted as dotted lines. A red-sequence model cor- 
responding z = 1.132 is represented as the solid black lines with a 
representative L* galaxy represented by the black diamond. 



luric line removal were performed using the well -exposed 
continua of spectrophotometric standard stars (Wade & 



Horne|[1988 [Foley et al.|[2003 |. 

i'hree independent redshitt determinations were per- 
formed using a cross-correlation al gorithm (IRAF 



RVSAO package; [Kurtz fc IVIink 1998[ ), a template fit- 
ting method (SDSS early-type F(JA templates), and a 
minimization technique by comparing to galaxy tem- 
plate spectra. There were only minor differences in the 
final results from the three methods. In total, we have 
obtained secure redshifts, consistent with membership in 
a single cluster, for 18 galaxies. Two of these galaxies 
have obvious [O II] emission, while the others have SEDs 
consistent with passive galaxies with no signs of ongoing 
star formation. 

A 3-(T clipping was applied around the peak in redshifts 
to select spectroscopic cluster members. Representative 
spectra of cluster members and a redshift histogram of 
cluster members are presented in Figure [3j Redshift in- 
formation for cluster members is presentedin Table [Tl A 
single galaxy was observed and has a secure redshift from 
both JVIagellan and VLT. Although the VLT spectrum 
shows clear Ca H&K absorption lines and the Magel- 
lan spectrum only shows the D4000 break, the measured 
redshifts are consistent. 

A robust biweight estimator was applied to the 
spectroscopic sample to determine a mean redshift of 
z = 1.131^q'oq3 and a velocity dispersion of — 

1230j^^gQ km s~^. The uncertainty in both quantities 
is determined through bootstrap resampling. Since the 
dynamics of passive and star-forming galaxies within 
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Fig. 3. — Left: Subset of optical spectra of galaxies identified to be cluster members redshifted by the cluster redshift, z = 1.132. The 
black vertical dashed lines show the position of [O II] A3727 and Ca H&K AA3933, 3968, while the grey dotted lines show the where the 
spectral features would fall at ±1000 km s"'^ relative to the cluster redshift. Right: Histogram of velocities for the spectroscopic cluster 
members relative to the biweight mean of all members. The blue and red histograms represent objects with and without obvious [O II] 
emission, respectively. 



clusters are different with passive galaxies assu med to 
be a better tracer of the cluster potential (e.g., |Fadda 



et al., .19961 


Girardi et al. 1996; 


ranyi & Ge 


ler 2000p, we also 



z = 1.132t°;g°^ and = 1270+^^;; km s-\ respectively. 

The velocity distribution of spectroscopic members rel- 
ative to the mean velocity (Figure |3| is skewed. The 
bootstrap resampling produced a positive skewness in 
89% of the samples, which indicates that it is not very 
statistically significant given our sample size. However, 
this may be an indication of substructure, consistent 
with th e po ssible merger seen in the X-ray image (see 
Section 2.4 1. There is no spatial segregation of galax- 



ies clustered in velocity space; however, the number of 
spectroscopic members is still small. 



2.4. X-ray Observations 

As part of a program to follow up the most massive, 
high redshift clusters detected with the SPT in 2009, 
we have obtained a 25 ks X-ray observation of SPT- 
CL J2106-5844. These data were obtained through a 
combined Chandra High Resolution Camera and Ad- 
vanced CCD Imaging Spectrometer Guaranteed Time 
Observation program. Data reduction was performe d 
with the method described by Andersson et al. (2010). 
The observation resulted in 1425 source counts within 
r5oo in the 0.5-7.0 keV band, corresponding to a flux, 
FxiO.5 - 2.0 keV) = (2.58 ±0.15) x lO'^^ erg cm-^ s^i. 
A 2" X 2" -binned X-ray image was produced using the 
X-ray events in the 0.7-4.0 keV band (Figure |4]). The 



TABLE 1 

Redshifts for. Galaxy Members of SPT-CL J2106-5844 



R.A. (°)'' 


Dec. (°) 




2^ 




Feature'^ 


Telescope 


316.46581 


-58.72418 


1 


1365 


(6) 


Ca H&K 


VLT 


316.47601 


-58.72635 


1 


1347 


(3) 


Ca H&K 


VLT 


316.48048 


-58.69972 


1 


1230 


(3) 


[O II] 


VLT 


316.48086 


-58.75017 


1 


1288 


(2) 


Ca H&K 


VLT 


316.49456 


-58.75092 


1 


1148 


(2) 


Ca H&K 


VLT 


316.50647<^ 


-58.73848 


1 


132 


(20) 


D4000 


Magellan 


316.50647'i 


-58.73848 


1 


1296 


(3) 


Ca H&K 


VLT 


316.51307 


-58.72857 


1 


1312 


(3) 
(6) 


Ca H&K 


VLT 


316.52563 


-58.75017 


1 


1450 


Ca H&K 


VLT 


316.53273 


-58.72584 


1 


120 


(20) 


D4000 


Magellan 


316.53346 


-58.77275 


1 


1335 


(5) 


Ca H&K 


VLT 


316.53434 


-58.75143 


1 


1389 


(4) 


Ca H&K 


VLT 


316.53679 


-58.75943 


1 


1274 


(6) 


Ca H&K 


VLT 


316.54011 


-58.71250 


1 


1384 


(3) 


Ca H&K 


VLT 


316.54123 


-58.76282 


1 


1384 


(5) 


Ca H&K 


VLT 


316.54297 


-58.76573 


1 


1260 


(4) 


Ca H&K 


Magellan 


316.54292 


-58.75800 


1 


1363 


(6) 


Ca H&K 


VLT 


316.55414 


-58.72558 


1 


1325 


(3) 


[O II] 


Magellan 


316.74972 


-58.72926 


1 


1136 


(5) 


Ca H&K 
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All coordinates arc in J2000. 
Uncertainties given in units of 0.0001. 

Spectroscopic feature primarily used in determining redshift. 
This is the same galaxy observed with both Magellan and VLT. 



observation reveals an X-ray luminous cluster with pro- 
nounced asymmetry in the surface brightness distribu- 
tion (see Figurepl, indicating that it may be in a merging 
state. In particular, we note that a bright concentration 
of gas is offset ~10"east of the main cluster centroid. 

X-ray observations of SPT-CL J2106-5844 were used 
to measure the gas temperature, Tx, and the X-ray ana- 
log of Comptonization, Yx = M^^s x Tx- The cluster 
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Fig. 4. — SPT-CL J2106-5844 at X-ray, millimeter, optical, and infrared wavelengths. Left: Chandra X-ray image from the 0.7-4.0 keV 
band, with 4x4 binningof the original ACIS pixels (0'.'49 on a side). The black contours correspond to the SZ significance from the 
left-hand panel of Figure [l| Right: LDSS3 optical and Spitzer/IRAC mid-infrared (;i[3.6] (corresponding to BRG channels) images. The 
frame subtends 4.'8 X 4.'8.T'he white contours are the smoothed X-ray signal from the left panel. 



mass can t hen be inferred with th e M-Yx and M-Tx 
relations of Vikhlinin et al. (2009). The rsoo radius is 
derived itera tively trom the M 



scaling relation as Santos et al. 



found that the cool-core fraction varies significantly be- 
tween different high redshift X-ray selected samples (e.g., 



described in Andersson et al. (2010). See Section |3.3 
for details regardmg the X-ray mass estimates. A spec 
trum is extracted from the data within the derived rsoo 
radius, excluding the core within 0.15 rsoo. The spec- 
trum is fit, keeping the absorbing hydrogen column fixed 
at the galactic value nn = 4.32 x 10^" cm~^, and the 
metal abundance fixed at 0.3 solar. Fixing the redshift 
to the mean spectroscopic redshift, z ~ 1.132, results in a 
best fit kT = n.Otl[l keV. Within rgoo (117"), including 
the cluster core, SPT-CL J2106-5844 has a luminosity of 
Lx(0.5 - 2.0 keV) = (13.9 ± 1.0) x lO"** erg s'^, compa- 



rable to the luminosity of the Bullet cluster (|Markevitch 
eFaL][2002| . 



An X-ray spectrum was produced from the data within 
20" (0.1 7 rsoo), where the metal abundance of the gas is 
likely to peak. The resulting spectrum is shown in Fig- 
ure[5]with a model folded through the detector response. 
The strong 6.7 keV Fe K line is present in the spec- 
trum at approximately 3 keV. This spectrum indicates 
z = 1.18 ± 0.03, a metal abundance Z = 1.44l[];^'j Zq, 

and a central temperature kT = 6.5j^!['^ keV. Notably, 
this measurement was the first robust spectroscopic red- 
shift measurement for this cluster. 
Merger simulations sugg est that high-redshift coo l-core 



clusters should be rare (Gottlober et al. 2001), and 
previous observations have mdicated a steep decline in 
the number of cool-co re clusters at high redshift (e.g., 
Vikhlinin et al. 2007). However, other analyses have 



2010). Future work comparing these re- 
sults to the cool-core fraction in SZ selected samples 
should help to understand the role of selection in these 
results. Regardless of the abundance at high redshift, 
SPT-CL J2106-5844 is one of the few high-redshift cool- 
core galaxy clusters known. 

3. MASS ESTIMATES FROM CLUSTER OBSERVABLES 

In estimating the mass of a galaxy cluster from any ob- 
servable property, care must be taken to properly account 
for Eddington bias, which, in this particular case, refers 
to the preferential scattering of lower-mass systems into a 
given observable bin due to a steeply falling source popu- 
lation. This is handled through application of Bayes' the- 
orem: we use scaling relations and noise measurements 
to calculate conditional probabilities of observable O at 
fixed mass, P{0\M), and combine these with a mass 
prior P{M) to yield a posterior estimate of the mass 
P{M\0). When the observable O is the same quantity 
that was used to detect the cluster, the prior probabil- 
ity of the cluster's mass is simply the cosmological mass 
function. However, when the observable in question is 
from a targeted follow-up observation, the prior must also 
take into account the selection function of the survey in 
which it was discovered. For follow-up observations tar- 
geting a single, very rare or interesting object, such as 
SPT-CL J2106-5844, the appropriate prior mass distri- 
bution is strongly influenced by the selection. For this 
reason, we report masses derived from follow-up obser- 
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12 5 12 

Energy (keV) Energy (keV) 

Fig. 5. — Core-excluded (circular annulus spanning 0.15-1 rsoo) and central region (circular aperture with radius 0.17 rgoo) X-ray spectra 
of SPT-CL J2106-5844 (left and right panels, respectively). Overplotted (solid line) is a fit APEC spectrum redshifted to 2 = 1.18 and 
folded through the detector sensitivity. The 6.7 keV Fe K line observed at ~3 keV is both obvious and strong for the spectrum from the 
central region. 



vations with a flat mass prior only. We report both flat- 
prior and mass- f unct ion-prior masses for the SZ-derived 
mass. In Section [3. 4[ we combine the conditional proba- 
bilities derived from different observables and apply the 
mass-function prior to the result in order to determine 
the best mass estimate based on all the available infor- 
mation. For the cosmological analyses in Section |4j we 
use only this joint posterior mass probability density. 

3.1. SZ Mass Estimate 
An SZ-based cluster mass est imate is produced fol- 



lowing [Vanderlinde et al.| ([2010| , using the single-band 
(150 GHz) detection significance ^ as a mass proxy. We 
direct the reader to that work for details. Briefly, an 
SZ-mass scaling r elation, derived from the simulations of 
Shaw et al. ( |2009 ) and constrained though a cosmological 
analysis, is used to estimate a set of conditional probabil- 
ities P(^|M), taking into account the measurement noise 
and the intrinsic scatter in the scaling relation. A poste- 
rior mass estimate is then calculated through application 
of a mass prior. The systematic errors are estimated by 
linearly expanding the mass estimate as a function of the 
scaling relation parameters and using the marginalized 
uncertainties on these parameters. 

The SZ detection of this cluster is from an untargeted 
survey, allow i ng us to take the halo mass function of 
Tinker et al. (2008) as a prior on its mass. We assume 



the WMAF 7- year preferred ACDM parameters (Lar 
son et al. 2010 ) when calculating this prior; marginaliz 
mg over the ACDM cosmological parameter uncertain- 
ties was found to have a negligible impact on the ef- 
fect of this prior. SPT-CL J2106-5844 has a single- 
band 150 GHz significance of ^ = 18.5, which results 
in an unbiased posterior SZ-derived mass estimate of 
Msz,2oo = (1.06 ±0.23) x lO^ ^ h^^ Mq 



We also report (see Section 3.4 and Table [2]) a biased 
SZ mass estimate, using a flat mass prior; we present it 
here for completeness only. It is important to note that 
this mass is not directly comparable to other flat-m ass- 
prio r est imates (e.g., the masses estimated in Sections 3.2 



and 3.3). The SZ observation was drawn from an untar- 
geted survey and hence from the entire underlying popu- 



lation of clusters (given by the halo mass function), while 
the targeted observations were conducted with a priori 
knowledge that a massive cluster was present; these two 
different modes of observation draw from different popu- 
lations and thus will exhibit different levels of Eddington 
bias. To reiterate, naively using any mass estimate with- 
out understanding the selection of the observed popula- 
tion can result in biased measurements. 

3.2. Dynamical Mass Estimate 

The dynamical mass of SPT-CL J2106-5844 can be 
estimated fro m the velocity disper sion using the M-ay 
calibration of Evrard et al. (2008). Numerical simula- 
tions indicate that galaxies are likely to exhibit veloc- 
ity bias compared to the dark matter v elocity field (e.g., 
Biviano et al.||2006[ [White et al.||2010[ ). Mock observa- 
tions of simulated clusters that model the galaxy selec- 
tion within followup observations like those reported here 
enable these biases to be quantified and corrected (Saro 
et al., in prep). Pubhshed results suggest these biases 
are a t the level of 10% in the dispersion (White et al. 
20101. Here we attempt no correction to the dispersion 
but account for this additional systematic uncertainty in 
the mass estimator. 

As the dynamical mass estimate arises from a targeted 
follow-up observation, the underlying mass distribution 
can be well-defined (and accounted for) only in the con- 
text of the SZ-based selection. Suc h a j oint estimate of 
mass is presented below in Section |3.4[ here we report 
an independent (hence biased) dynamical mass estimate 
using a flat prior on the mass. 

Using the velocity dispersion from all 18 spectroscopic 

cluster members gives Mdyn,200 — l-Sloi ^ ^^^^ ^70 ho- 
using the velocity dispersion from only the passive galax- 
ies provides a very similar result, Mdyn,200 = l-4lo;8 ^ 
1015 /j-i Mq). In both cases, the uncertainty takes into 
account the statistical uncertainty of the velocity disper- 
sion, an additional 13% intrinsic scatter of the velocity 
dispersion due to the orien tation of the veloc ity ellip- 
soid along the hue of sight (|White et al.||2010[ ), the 4% 

( |:jO(m^ -a„ rela- 



intrinsic scatter of the ^Evrard et al 
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tion, and a 10% estimate of the systematic uncertainty 
in the dispersion that arises from the currently uncer- 
tain kinematic relationship between the galaxies and the 
underlying dark matter. 

3.3. X-Ray Mass Estimate 
As described in Section |2.4| the X-ray observations 



of SPT-CL J2106-5844 were used to measure the gas 
temperature, Tx, and the X-ray analog of Comptoniza- 



tion, Yx 



X T: 



Following the methodology 
of Andersson et al. (20101, the data were iteratively 
fit to deter mine r..;nni '^'y/and Yx - Using the scaling 



(2009), these values pro- 



relations of Vikhlinin et al. 

vide estimates ol the gas and total mass of the cluster, 
Mgas,500 = (1.27 ± 0.07) X 1014 h^^ M© and My^^goo = 
(0.93±0.19) X 10^^ hjQ Mq, respectively. Estimating the 
cluster mass from the gas temperature, we find a simi- 
larly high mass of Mtx, 500 = (0.92 ±0.37) x 10^^ /ifq M©, 
where we ha ve increased the mass by 17% according to 
the recipe in Vikhlinin et al. ( |2009 ) because of the large 
degree of asymmetry m the surface brightness distribu- 
tion, indicative of merging activity. 

For both mass estimates above we have adopted a sys- 
tematic uncertainty of 9% for the mass calibration of 
the scaling relations. The uncertainty on the self-similar 
evolution of these relations was estimated using the sim- 
ulation informed estimate of 5% and 7% uncertainty on 
the normalization at z = 0.6 for the M-Yx and M-Tx, 
respectively. This uncertainty on the evolutionary term 
was extrapolated to the cluster re dshift. For further dis - 
cussion on these uncertainties, see Vikhlinin et al. ( 2009 ); 
the appropriate mass prior to account tor Fiddington bias 
is only well-defined in the context of SZ-based selection. 

3.4. Combined Mass Estimate 

The X-ray, SZ, and optical redshift data can each 
be used to produce an estimate of the mass of SPT- 
CL J2106-5844. We list our estimates for the mass of 
SPT-CL J2106-5844 from each of these observables in Ta- 
ble [2] Note that we have converted the My-^^-^soo mass to 
Myx,200 assuming a Navarro-Frenk- White densi ty pro- 
file an d the mass-concentration relation of , Duffy et al 
(2008). Mtx, 500 is not used in the combined mass esti- 
mate since it is correlated with My^ ,500 and has a larger 
uncertainty. Similarly, Mdyn,200 is not used in the com- 
bined mass estimate because of its large uncertainty. 

We assume for the purposes of the combined mass es- 
timate that the uncertainties are uncorrelated between 
the two masses. In principle, the different observables 
could have correlated scatter in their scaling with mass, 
however, this is expected to be negligible at the current 
S /N of the observables. This simplification allows us to 
write the posterior mass distribution as a product of the 
two independent conditional probabilities. 



P{M\C,Yx) (X P{M) ■ Pirn) ■ P{Yx\M), (1) 



where P(M) is again the halo mass function of Tin- 
ker et al.| ( |2008[ ), evaluated at the WMAP 7-year pre- 
terred point m ACDM parameter space. Functionally, 
this reduces to a product of two probability density 
functions, the unbiased SZ mass estimate and the fiat- 
prior Yx mass estimate. We find the combined, un- 



TABLE 2 

Mass Estimates for SPT-CL J2106-5844 



Observable 



Measurement 



Af200(10l5?t-iMQ) 



SZ ^ 18.5 
SZ 5 (fiat prior) 18.5 



1.06 ±0.23 

1.24 ±0.30 



Yx (14.3 ± 3.4) X IQi'' Mq keV 1.85 ± 0.40 

Tx ll.OtJ J keV 1.83 ±0.76 

o-t,^ 1270^220 ^™ '^"^ ^■'^-ol 

Combined ■■■ 1.27 ±0.21 

Note. — Unbolded masses indicate Eddington-biased mass estimates, 
calculated using flat priors on mass. Note that the SZ (untargetcd) 
measurement suffers from a considerably different Eddington bias than 
the other (targeted) flat-prior estimates. The probability distribution 
for the combined mass was obtained by multiplying the unbiased (bold) 
SZ mass probability by the flat-prior Yx mass probability (see text for 
details). Only the unbiased SZ ^ and Yx mass estimates were used to 
generate the combined mass estimate. 

^ As determined by the biweight estimate of the passive galaxies. 

biased, mass estimate from the two observables to be 
Afaoo = (1.27 ± 0.21) x lO^^ hj^^ M©. 

All mass estimates are consistent and very high for a 
cluster at this redshift. Each individual mass estimate, 
as well as the combined estimate, indicates that SPT- 
CL J2 106-5844 is the most massive cluster yet discovered 
at z > 1. It is ^60% more massive than SPT-CL J0546- 
5345 ( [Brodwin et al.||2010[ ) at z = 1.07, the second most 
massive galaxy cluster at z > 1. 



4. LIKELIHOOD ANALYSIS 

Since SPT-CL J2106-5844 is the highest S/N clus- 
ter found in 2008-2009 SPT survey fields, the highest- 
redshift spectroscopically confirmed SPT cluster, and the 
most massive cluster discovered at z > 1 by any method, 
it raises the question: Is the existence of SPT-CL J2106- 
5844 consistent with a Gaussian-distributed primordial 
density field in a ACDM universe? 

We investigate the probability of finding a cluster with 
the mass and redshift of SPT-CL J2106-5844 in the con- 
text of AC DM as follows. Based on the likelihood cal- 
culation of 'Vanderlinde et al. (2010) and the CosmoMC 
software package ( jLewis fe Bridle ^2002) , we produce an 
MCMC chain sampling the parameter space of spatially 
flat ACDM cosmolo gy, storing the halo mass function of 
Tinker et al. (2008) at each step in the chain. Current 
data sets, c onsisting of the WMAP 7-year data ('Larson 



et al.||2010|) and t he SPT 2008 cluster catalog dVander 



linde et al. 2010), are used to constrain the parameter 



space explored. I'his yields a set of halo mass functions 
representative of current uncertainties of ACDM param- 
eter space. 

Integrating each of these mass functions, we compute 
the Poisson likelihood of finding a cluster at or above the 
mass and redshift of SPT-CL J2106-5844 in a 2500 deg^ 
field. For the MCMC chain described above, the ex- 
istence of a cluster at or above the mass and redshift 
of SPT-CL J2106-5844 is less than 5% likely in mod- 
els contributing 32% of the WMAP7+V10 probability. 
Marginalized over all WMAP7+V10 allowed parameter 
space, we find that there is a 7% chance of finding a clus- 
ter at least as massive and at a redshift at least as high 
as SPT-CL J2106-5844 in the SPT survey field. Prior to 
having observed this cluster, the expectation would be 
that there is only one galaxy cluster at least as massive 
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and at a redshift at least as high as SPT-CL J2106-5844 
in the entire sky. 

5. CONCLUSIONS 

We present the discovery and initial characterization 
of SPT-CL J2106-5844. This massive cluster was discov- 
ered in the 2500 deg^ SPT survey as a highly significant 
SZ decrement. Initial follow-up in the optical and in- 
frared indicated that it was both at high redshift and 
very rich. Optical spectroscopy yielded a cluster red- 
shift of z — 1.132, consistent with the redshift found 
from the X-ray spectrum of the hot intracluster gas, 
z = 1.18, and the broad-band SEDs of member galax- 
ies. SPT-CL J2106-5844 is currently the highest-redshift 
SZ-discovered galaxy cluster. 

SPT-CL J2106-5844 is also extremely massive. X- 
ray, SZ, and velocity dispersion data provide comple- 
mentary and consistent measurements of the mass (see 
Table [2]). The best combined estimate of the mass is 
M200 = (1.27 ± 0.21) X 10^5 h^^ Mq, making SPT- 
CL J2106-5844 the most massive galaxy cluster yet dis- 
covered at z > 1. The combination of its high mass and 
high redshift make SPT-CL J2106-5844 a very interesting 
object in which to study galaxy formation and evolution 
only a few Gyr after the formation of the first stars. The 
X-ray data show that SPT-CL J2106-5844 has a cool 
core and may be currently undergoing a merger. The 
distribution of galaxy velocities is slightly skewed, but 
is consistent with a Gaussian distribution. Additional 
spectroscopy may help determine if there is significant 
substructure in the cluster. 

The high mass and redshift of SPT-CL J2106-5844 also 
make it an interesting object from a cosmological per- 
spective. The chances of finding a cluster this extreme 
in mass and redshift in a 2500 deg^ survey is only 7%, 
given current cosmological parameter constraints, with 
only one such cluster expected in the entire sky. As the 
SPT cluster catalog grows, and our understanding of the 
various mass-observable scaling relations improves, SZ- 
discovered clusters will enable yet more powerful tests of 
the current cosmological paradigm. 

Facilities: Blanco (NEWFIRM), CXO (ACIS), Mag- 
ellan:Baade (IMACS), Magellan:Clay (LDSS3), Spitzer 
(IRAC), South Pole Telescope, VLT:Antu (F0RS2) 
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